TIFR/TH/05-40 



Dynamics of Fluctuation Dominated Phase Ordering: Hard-core Passive Sliders on a 
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We study the dynamics of a system of hard-core particles sliding downwards on a one dimensional 
fluctuating interface, which in a special case can be mapped to the problem of a passive scalar 
advected by a Burgers fluid. Driven by the surface fluctuations, the particles show a tendency to 
cluster, but the hard-core interaction prevents collapse. We use numerical simulations to measure 
the auto-correlation function in steady state and in the aging regime, and space-time correlation 
functions in steady state. We have also calculated these quantities analytically in a related surface 
model. The steady state auto-correlation is a scaling function of t/L z , where L is the system size 
and z the dynamic exponent. Starting from a finite intercept, the scaling function decays with a 
cusp, in the small argument limit. The finite value of the intercept indicates the existence of long 
range order in the system. The space-time correlation, which is a function of r/L and t/L z , is 
non-monotonic in t for fixed r. The aging auto-correlation is a scaling function of ti and ti where 
t\ is the waiting time and t% the time difference. This scaling function decays as a power law for 
ti ti; for ti S> ta, it decays with a cusp as in steady state. To reconcile the occurrence of strong 
fluctuations in the steady state with the fact of an ordered state, we measured the distribution 
function of the length of the largest cluster. This shows that fluctuations never destroy ordering, 
but rather the system meanders from one ordered configuration to another on a relatively rapid 
time scale. 



I. INTRODUCTION 

The concentration of a substance advected by a driving 
field such as a fluid flow often shows interesting behavior. 
In examples such as smoke dispersed in air or fluorescent 
dye carried by a turbulent jet, an initial local concen- 
tration of the passive particles generally spreads out in 
space 0. However, if the fluid is compressible, instead of 
spreading out, the advected substance may show a ten- 
dency to cluster, as with air bubbles in water or dust 
particles in air 0, . 

In the situations discussed above, the advected sub- 
stance has a negligible effect on the fluid flow. These 
problems are examples of 'passive scalar' problems, where 
the dynamics of a non-equilibrium driving field strongly 
affects that of the other (passive) field with no back-effect 
from the latter. In this paper, we study clustering in a 
set of passive particles which are subjected to a fluctuat- 
ing force field. The specific system we study consists of 
hard-core particles sliding under gravity on a one dimen- 
sional fluctuating interface; the instantaneous force on a 
particle is then proportional to the local slope of the sur- 
face. Interestingly, our results also pertain to particles 
advected by a fluid, using the fact that the equation gov- 
erning a moving interface can be mapped onto the Burg- 
ers equation, which describes a compressible fluid 4'j. 

The degree of clustering of the particles depends 
strongly on the interactions between them. Two cases 
have been studied earlier — particles which are completely 
non-interacting 0, 01 and particles which interact via 
hard-core exclusion [3, @. We study the latter, more 
realistic, case in this paper and address two broad issues, 
namely dynamics and ordering. 

The first issue concerns the time-dependent proper- 
ties of the passive particles. We obtain results for the 



auto-correlation and space-time correlation functions in 
steady state, and aging correlations in the approach to 
steady state. To our knowledge, the dynamics of pas- 
sive scalars has not been explored systematically, and 
our study adds to the relatively sparse work on this im- 
portant question We have considered two kinds of 
surface evolutions — namely those that respect symme- 
try under reflection (Edwards- Wilkinson type) and those 
which break this symmetry (Kardar-Parisi-Zhang type). 
Our simulation results for the dynamics of the sliding par- 
ticle system are supplemented by analytical calculations 
for a related model of coarse-grained interface variables. 



The second issue concerns the characterization of the 
steady state as an ordered state. Earlier studies of static 
properties revealed that strong fluctuations are present 
in the steady state and do not decrease even in the ther- 
modynamic limit 0. On the other hand, the scaling 
function for the density correlation indicates that the 
system has long range order. This sort of fluctuation 
dominated phase ordering (FDPO) is characterized by 
a broad distribution of the order parameter. The ques- 
tion then arises: In what sense does FDPO represent 
an ordered state, if strong fluctuations drive it between 
macroscopically different configurations on a relatively 
rapid time-scale? We address this by studying the varia- 
tion of the length of the largest particle cluster present in 
the system and show that the corresponding probability 
distribution provides an unequivocal signal of ordering. 
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II. OVERVIEW 
A. Surface Fluctuation and Particle Movement 

A surface with no overhangs is completely specified by 
the height h(x,t) at point x at time t. The evolution of 
the height field is taken to be described by the Kardar- 
Parisi- Zhang (KPZ) equation Q. 

dh d 2 h , fdh\ 2 , . 

The first term represents the smoothening effect of sur- 
face tension v\, and r)i(x,t) is a white noise with zero 
average and (r)i(x, t)rjx(x', £')) = T5(x — x')5(t — t'). No- 
tice that if A ^ 0, h ^ —h symmetry is not preserved, 
reflecting the fact that the interface moves in a preferred 
direction. However, if A = 0, the equation has an h — ► —h 
symmetry and describes the Edwards- Wilkinson (EW) 
model Hg. 

The height-height correlation function has a scaling 
form for large separations in space and time |ll| : 

([h(x,t) - h{x',t')f) « \x-x'\ 2 *f (j^^jj) (2) 

Here / is a scaling function and \ and z are the roughness 
and dynamic exponents, respectively, with values which 
depend on the surface dynamics. For an EW interface 
X = 1/2, z = 2 while for a KPZ interface x = V 2 , 
z = 3/2. 

The hard-core particles slide downwards along the lo- 
cal slope (tj!) of the interface. In the overdamped limit, 
their velocity is proportional to the local gradient of 
height. The equation governing the evolution of parti- 
cle density can be derived from the continuity equation 
dp( ^ t ,t ^ = — d* 7 ^/*" 1 . The local current J(x,t) has a sys- 
tematic part p(l — p)(l — 2t^), a diffusive part — ^2§§ 
(driven by local density inhomogeneity) and a stochas- 
tic part rj 2 (x,t) (a Gaussian white noise). The time- 
evolution equation for the density fluctuation p = p — po 
is then 

f ^ 2 g + 2p o (l-Po)0 
-(l-2p -2p)(§) [l-2(f)] + 
2(l-2p o )p0-2p30 + ^ (3) 

We will not analyze this equation directly; rather we 
will study the particle dynamics by performing numer- 
ical simulations on a lattice model, described in section 
3, whose long distance and long time properties are ex- 
pected to be described by Eqs.J2J) and ©. 

B. FDPO : Static Properties 

In an earlier study on static properties of this model Q , 
the density-density correlation C (r, L) of the sliding par- 



ticles was measured and found to be a scaling function of 
r/L in the scaling limit r — > oo, L — > oo with r/L fixed, 
as for phase-ordered states. The scaling function has a fi- 
nite intercept m 2 , and for small argument it decays with 
a cusp: 

C(r,L) = /(I) (4) 

m 2 is a measure of long range order (LRO) as LRO is 
defined by the large r (scaling limit) behavior of the cor- 
relation function. In the limit of an infinite system, this 
corresponds to r/L — > 0. 

The value of the cusp exponent a depends on the dy- 
namics of the driving surface: a ~ 0.5 for particles on an 
EW interface while a ~ 0.25 when the underlying surface 
is of KPZ type 0. In a related coarse-grained surface 
model defined in Q and discussed in section 3 below, the 
correlation function was shown analytically to have the 
above scaling form with m 2 = 1 and a = 0.5 for both EW 
and KPZ surfaces. For hard-core particles sliding on a 
two dimensional surface, the same scaling form is found 
for C(r,L), but with a different value of the intercept 
and the cusp exponent. For customary phase-ordering 
systems, it is expected that this scaling function should 
decay linearly, consistent with the Porod law The 
cusp (a < 1) is one manifestation of the unusual nature 
of FDPO. 

For a phase-ordered system, the lowest non-zero 
Fourier component of the density profile measures the 
extent of phase-separation and is an appropriate order 
parameter. A characteristic of FDPO is that the distri- 
bution of this order parameter remains broad even in the 
thermodynamic limit, indicating the presence of strong 
fluctuations. 



C. FDPO : Dynamical Properties 

In this paper, we study the dynamics associated with 
FDPO. Our results are summarized below. We find 
that the steady state auto-correlation in the density fluc- 
tuation a(x,t) of the sliding particles (SP): Asp(t) = 
(a(x,0)a(x 1 t)} is a scaling function of t/L z , consistent 
with the notion of phase ordering. For t <C L z , this 
scaling function decays with a cusp, in contrast to a lin- 
ear decay normally expected for phase-ordering systems. 
The presence of the cusp is the dynamical manifestation 
of the unconventional character of FDPO. 

We also monitored the auto-correlation function in the 
aging regime Asp(h,t 2 ) = (a(x,ti)a(x,ti +t 2 )). From 
the dynamic scaling hypothesis, Asp{ti,t 2 ) should be a 
scaling function of t\ /t 2 , the ratio of initial time to the 
time lag [l^. In the limit t 2 ti, the scaling function 
shows a power law decay with an exponent that depends 
on the phase-ordering kinetics. In the opposite limit t\ 3> 
t 2 , the scaling function decays with a cusp as in steady 
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state, the only difference being that the system size L, as 
appears in the steady state scaling function, is replaced 

by the coarsening length t^ z . 

The space-time correlation for the sliding particles in 
steady state, defined as Gsp{r,t) = (a(x, 0)a(x + r, t)) is 
a function oir/L and t/L z . When plotted against t/L z 
for a fixed value of r/L, this function is non-monotonic. 
and decays like the steady state auto-correlation for large 
t. The auto-correlation function in steady state and in 
the aging regime, together with the space-time correla- 
tion function in steady state constitute our dynamical 
characterization of FDPO. 



D. FDPO : An Ordered State? 

The distinctive feature of FDPO is the presence of 
strong fluctuations in steady state. On the one hand, a 
diagnostic of LRO such as a non-zero value of the inter- 
cept m 2 indicates an ordered state. On the other hand, 
the macroscopic state of the system changes relatively 
rapidly over a time-scale ~ L z , in contrast to a normal 
phase-ordered system, where the typical time grows ex- 
ponentially with L. The small lifetime of a given macro- 
scopic ordered state gives pause, and seems to contradict 
the notion of LRO. To resolve this apparent contradic- 
tion we study the time dependence of the length of the 
largest particle cluster l max (t) . We conclude that the 
short lifetime is associated with the system wandering 
over a multitude of ordered states, each very different 
from the other, but all characterized by a large value of 
lmax(t). Dynamical excursions away from this attractor 
of ordered states are extremely infrequent and associated 
with an exponentially growing time scale. 



III. DESCRIPTION OF THE MODEL 

We study a discrete model of a fluctuating interface on 
which hard-core particles slide downwards under gravity, 
following the local slope of the interface. The 1-d inter- 
face of length L, consists of discrete surface elements; the 
slope of the surface elements between the i-th. and + 
th site is r i+ i, which can take the value +1 or — 1. Ac- 
cordingly the height at site i is given by hi = J2]=i T j-\- 
The dynamics of the interface involves stochastic corner 
flips with exchange of adjacent r's; the transition /\ to \/ 
occurs with a rate p\ while \/ to /\ with rate q\. We have 
taken p\ = qi = 1 for an EW surface and pi = 1, qi = 
for KPZ surface. The overall slope T = Ylt=i T i+± 1S 
conserved and in our case we will consider T = 0, mean- 
ing that the interface is untilted. 

The hard-core particles are represented by variables 
{<jj} each of which takes a value +1 or —1 according 
as the i-th site contains a particle or a hole. The devi- 
ation from half-filling S = J^Li <J% is conserved and 
we consider 5 = 0, corresponding to the case of half- 



filling. A particle and hole on adjacent sites + 1) 
exchange with rates that depend on the intervening local 
slope T i+ i ; thus the moves »\o — > o\» and o/» — > »/o oc- 
curs at rate P2 while the inverse moves occur at rate qi. In 
the case when the particles are sliding downwards along 
gravity, we have qi < pi- We have considered p% = 1 
and qi = 0. Because of the hard-core exclusion between 
the particles, with the above update rules, the system 
has a particle-hole symmetry, i.e. any correlation func- 
tion involving the density variable remains invariant as 
the particle density is replaced by the hole density. This 
implies that unlike the case of non-interacting particles 
where the correlation function of the density variable is 
qualitatively different between the advection case and the 
anti-advection case 0, the two processes show identical 
correlation behavior here. 

From the dynamical rules, it follows that the movement 
of particles depends on the fluctuations of the underly- 
ing interface. Due to gravity the particles tend to slide 
down into local valleys. However, in the non-equilibrium 
system under consideration, before the particles can fill 
in the lowest valley, the interface evolves, often causing 
the valley to turn over. Nevertheless, it is useful to con- 
sider the adiabatic limit where the interface moves in- 
finitely more slowly than the particles, in which case the 
particles have ample time to explore the landscape and 
eventually settle in the deepest valleys. It seems plausi- 
ble that the dynamics of hills and valleys of the interface 
may provide insight into the dynamics of the particles. 
This motivates the definition of a coarse-grained depth 
model (CD model) as follows \t\. Consider the variable 
Si{t) defined as Sj(i) = sgn[hi(t) — (h(t))], where (h(t)) 
is the average height at time t: (h(t)) = 4 Ya=i hiit)- 
The variable Si(t) can take values +1,-1 or 0, depending 
on whether the position of the ?-th site is above, be- 
low or at the average level. In other words, Si(t) gives 
a coarse-grained description of the surface by labeling 
'highlands' and 'lowlands'. For an EW interface, the dy- 
namics is tractable and we obtain an analytic expression 
for time-dependent correlations of Sj(t). These results 
might be expected to be close to those of Ui{t) in the ex- 
treme adiabatic limit. As a matter of fact, we find that 
they also describe qualitatively the particle model even 
in the strongly non-equilibrium case. 



IV. AUTO-CORRELATION FUNCTION IN 
STEADY STATE 

We have studied the auto-correlation A(t, L) involving 
the density variable (e>i(0)o"j(t)) in the sliding particle 
(SP) model and also that for the CD variables (si(O)sj(t)) 
in the CD model. Periodic boundary conditions are used. 
We will see below that in the steady state of a system of 
size L, the auto-correlation A(t, L) is a scaling function 
of 4j, where z is the dynamic exponent defined earlier. 
This scaling function shows a cusp in the small argument 
limit, as seen previously in the static correlation scaling 
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CD Model 


SP Model 


EW 


KPZ 


EW 


KPZ 


m 


1.0 


1.0 


0.82 ±0.03 


0.75 ± 0.04 


P' 


0.25 


0.31 ±0.02 


0.22 ±0.02 


0.18 ±0.01 


A 


0.5 


0.5 


0.26 ± 0.005 


0.15 ±0.005 


7 


0.75 


0.84 ±0.02 


0.69 ±0.02 


0.82 ± 0.04 



TABLE I: The values of relevant exponents and order parame- 
ter for dynamical characterization of CD model and SP model 



function [Eq. ©] 



points goes as L~ x 14]. For small t, the probability that 
any one of these points will actually take part in a lo- 
cal fluctuation is proportional to t. This immediately 
implies Aco(t < 1, i) ~ 1 — b\j^. Comparison with 
Eq. © shows that for the CD model, we have A = \ = | . 
Note that although EW and KPZ surfaces have different 
dynamics, the above argument holds for both of them as 
their stationary measure is the same in l-d. 



0.98 



- w 



KPZ 



A(t,L) 



1 - b 



t 



(6) 
(7) 



m is a measure of the LRO as explained in section 2. Note 
that for large enough time, the auto-correlation function 
is expected to factorize and become equal to m 2 . 

However, for small time, t < 1, which falls outside 
the scaling regime, the auto-correlation function shows a 
linear drop with an L-dependent slope: 



A(t,L) 



(t<l) 



(8) 



If to 2 = 1, as shown below for the CD model, matching 
Eqs.© and © for t ~ 1 yields 



A = zf3' 



(9) 



If m 2 1, as happens for the SP model, a relation be- 
tween the exponents cannot be obtained. Instead, the 
matching condition determines a time scale t* for the 
crossover from the linear decay in Eq.© to the cuspy 
decay in Eq.Q. In the large L limit, we find to the 
leading order, 



t* = 



1 — m l 



-L 



(10) 



Let us illustrate these properties, by discussing the 
auto-correlation in the CD model, defined as AcD(t, L) = 
(si(O)si(t)}. First consider short times t < 1. At t = 
let the initial configuration of the surface be {hi(0)}. 
As time passes, there are stochastic corner flips, as de- 
scribed in section 3. However, only those flips occur- 
ring close to the average level can cause a change in the 
CD variable Sj(i), as any local fluctuation far above or 
below the average level, would not change the sign of 
Si(t) = (hi(t) — (h(t))). More precisely, only those sites 
in {hi(0)} which have at least one neighbor situated ex- 
actly on the average level, putatively contribute to the 
drop in auto-correlation function. Now, for a self-affinc 
surface of length L and roughness exponent \i the num- 
ber of such points scales as L 1_x and the density of such 



0.92 



EW 



0.86 



0.2 



0.4 



t/L 



FIG. 1: Illustrating the linear drop of A(t,L) for short times 
t < 1 in the SP model for system size L — 128, 256, 512. 

For the particle model, although the initial drop is 
found to be linear as described in Eq.©, the exponent 
A takes the value 0.26 ± 0.005 for particles on an EW 
surface and 0.15 ± 0.005 for particles on a KPZ surface. 
The data are shown in Fig.©. 

For t > 1, we have analytically calculated Aco(t,L) 
for an EW interface. This exploits the fact that hi(t) in 
this case is a Gaussian variable, implying Sj correlations 
satisfy the arc-sine law 



(sj(t)sj(O)) = - sin" 

7T 



(Hi(t)Hi(0)) 



(11) 



where Hi(t) = hi(t) — (h(t)), which is also a Gaussian 
variable. If hk(t) is the Fourier transform of hi(t), the 
numerator in the argument of arcsine can be written as 

Ek^(hk(t)h-k(0)) = E^l 3 ^, using the dis- 
crete version of the EW equation. Here, Ck = Av\ sin |. 



Moreover, {H.?(t)) 
have 



k^a-k- Thus we 



(si{t) Si (0)) = -sin 

TT 



E cxp(-c t t) 



(12) 



We have numerically evaluated this discrete sum and 
plotted it in Fig. (2a) against the scaling argument t/L 2 
for different L values. The cusp exponent can be read off 
from the plot in the inset. 
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In the continuum limit, Eg. 1)12(1 becomes 



(s(x,t)s(x,0)) = — sin 

7T 



fldk 



— p — 



n dk 



(13) 



The 

LF 



integral in the 

(-^)+V^er/ 



numerator 

27rVf 
L 



takes the form 



In the 



limit t/L 2 <C 1, this becomes, to the leading order, 
IT _ VTrf]. Noting that the denominator is ^--^ 
and expanding for small values of we get 

1/4 



(s(a:,t)«(a;,0)) 



4 / t 



(</i 2 < 1) (14) 



Comparing with Eq.Q gives m 2 — 1, /?' = i, z = 2. 

For the KPZ surface, the time evolution equation for 
the height field is not Gaussian and hence such an ana- 
lytical treatment is not possible. We study AcD(t) using 
Monte Carlo simulation. No initial equilibration is re- 
quired as the steady state measure for a KPZ surface with 
periodic boundary conditions gives equal weight to every 
configuration. We followed the update rules discussed in 
section 2 and averaged over sites as well as over 10 5 his- 
tories. The results are shown in Fig. (2b). A good scaling 
collapse is obtained for different L, on rescaling the time 
to t/L z with z = §. The cusp exponent (3 1 was extracted 
by plotting m 2 — Aco{t) against t/L z (shown in the in- 
set), using m 2 — 1 and this gives (3' to be 0.31 ± 0.02. 
Our best estimate corresponds to the largest system size 
L = 2048. The error-bar is based on the values of (3' 
obtained for smaller system size (L = 512,1024); the 
statistical error is much smaller. 

For the sliding particle (SP) model, the steady state 
measure is not known analytically. In our simulation, we 
started from a randomly disordered configuration and al- 
lowed a long time ~ 10L Z for the system to reach a steady 

state. We then measured J2i=i cr i(0) cr i(^) f° r approxi- 
mately L z time-steps. We waited several thousand time- 
steps before repeating the procedure, and averaged over 
10 4 histories . 

For particles sliding on an EW interface, we obtained 
a good scaling collapse of A$p(t, L) for different L after 
rescaling time to t/L 2 [Fig. (2a)]. The cusp exponent was 
extracted by fitting m 2 — Asp(t, L) to a power law. We 
have estimated m 2 by using the same technique as dis- 
cussed in |7J . The best estimate of m 2 corresponds to the 
value for which the structure factor has the largest power 
law stretch. We found that m 2 shows a systematic de- 
pendence on L and the cusp exponent j3' is in fact quite 
sensitive to the value of m 2 . We have used m 2 ~ 0.82, 
our estimate from the largest system size we could ac- 
cess (L = 4096). This yields /3' ~ 0.22. On the other 
hand, using m 2 ^ ~ 0.85, which we get by extrapolating 
the dependence of m 2 on L for an infinite system, we find 
(3' ~ 0.20. 

For the SP model on a KPZ surface, we find m 2 ~ 
0.75. Figure (2b) shows the scaling collapse for different 




FIG. 2: Scaled auto- correlation function in steady state for 
SP and CD models for (a) EW and (b) KPZ interfaces. In 
both cases, we used L=512, 1024,2048. The cusp exponents 
were estimated using the plots shown in the inset. 



L after rescaling the time by I? 1 2 . The inset shows that 
m 2 — Agp(t, L) follows a power law and the exponent is 
found to be j3' ~ 0.18. The value of f3' obtained using 
m 2 is ~ 0.17. 



V. AUTO-CORRELATION IN AGING REGIME 

The aging auto-correlation function A(ti, £2) is defined 
as {ai{t\)ai{t\ -M2)) for the particles and as {si(t\)si{t\ + 
t 2 )) for the CD variables. A(ti,t 2 ) depends on both t\ 
and t 2 . For 1 <C t 1 ,t 2 <C L z , ^(^^2) is a function of 
1^-, as expected for phase ordering systems [l3- I n the 
limit when ti 3> t%, this scaling function has a power law 
decay 



(15) 



-4(*i,fe)~ ( ^ ) fori 2 >ii, 



while in the opposite limit, t\ S> t 2 , the scaling function 
has the form 
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This is similar to the form of the steady-state auto- 
correlation in Eq.(j2J) with L replaced by t\ , meaning 
that locally the system has reached steady state over a 
length scale of t\' z . 



We first present our results on the CD model. As in the 
case of steady-state auto-correlation, we have been able 
to calculate Aco{ti,t2) for an EW surface analytically. 
Following similar steps as in the last section, we obtain 



AcD(ti,t 2 ) = —sin' 

7T 



E 



CXp( — Cfct2)-CXp[-Cfc(2tl+t2)] 

Cfc 



r v i-c XP (-2c fc ,t 1 ) i 1/2 ry i 



-cxp[-2c fc „(t 1 +t 2 )] 



1/2 



(17) 



Taking the continuum limit and using t\,t 2 <C L 2 , we 
obtain 



VI. SPACE-TIME CORRELATION IN STEADY 
STATE 



AcD(ti,t 2 ) = -sin' 
ir 



y/2t! + t 2 



(2t 1 ) 1/i (2t 1 + 2t 2 ) 



1/4 



(18) 



3/4 



In the limit t 2 ^> t±, right hand side becomes ^ {^) 

Comparing with Eq. l|15[l . we get 7 = §. In the opposite 
limit, when t% 3> t 2 , the right hand side becomes, after 
simplification, 



A CD (ti,t 2 ) « 1 



2t 

7T 



1/4 



(19) 



Comparing with Ea. (|16fl . we find (3' — 1/4, as expected. 

Figure (3a) shows the numerical evaluation of the dis- 
crete sum in Ea. H17|) . The power law characterizing the 
decay has been shown in the inset. 

In our Monte Carlo simulations, we have a spatial av- 
erage as well as an average over 10 4 histories. For the 
CD model of a KPZ surface, we started with a flat in- 
terface as an initial condition and evolved it in time to 
measure AcD{ti, t 2 ). The results are shown in fig. (3b). 
The best estimate of the cusp exponent corresponds to 
t\ = 32000 and the error bar is based on its values for 
h = 2000,8000. This finally gives f3' = 0.31 ± 0.01, 
which is close to the steady state value. The inset shows 
the power law decay and the exponent 7 takes the value 
0.84 ± 0.03. Here, the best estimate is for t\ = 500 and 
the error-bar is for t x = 2000, 8000. 

For the SP model on an EW interface, we start with 
randomly distributed particles on a random surface pro- 
file. The aging auto-correlation Asp(ti, t 2 ) shows a scal- 
ing collapse as plotted against £2/^1 [see fig. (3a)]. The 
value of the cusp exponent (3' is 0.20 ± 0.02, close to 
its steady state value. The inset shows plot in the 
regime t 2 3> t\. The power law exponent in this case 
is 7 = 0.69 ±0.02. 

The SP model on a KPZ surface also starts with the 
random initial condition. The data are shown in fig. (3b). 
The exponents are j3' = 0.17 ± 0.01 and 7 = 0.82 ± 0.04. 



In this section, we discuss the behavior of space- 
time correlation G(r, t, L) defined in steady state as 
(ai(0)<T i+r (t)) for the particles and (si(0)si +r (t)) for 
the CD variable. G(r, t, L) does not show any 
L— independent scaling between r and t. Rather, it is 
a function of the scaled variables £ = r/L and t = t/L z 



G(r,t,L)=g({,T). 



(20) 



With £ held fixed, g shows an interesting non-monotonic 
behavior with r. <?(£, 0) reduces to the pair correlation 
function /(£) [see Eq.(@J], and as r increases, g(£,,r) is 
observed to rise and attain a peak [see Fig. (4a)]. Fi- 
nally for larger r, it decreases and merges with the auto- 
correlation scaling function h{r) [see Eq.(|HJl]. Note that 
G(r = 0,t,L) = A(t,L), and by continuity, we ex- 
pect that when £ z <§; r, the scaling function should be- 
have like h(r). From our knowledge of the scaling func- 
tions /(£) and h(r), we have been able to verify that 
/(£) < h(r = £ z ). This implies that g(£, r) must show 
an initial rise. 

For the CD model on an EW interface, 



GcD(r,t,L) = —sin 1 

7T 



E 



fc>0 



cxp( — Cfct)2 cos(kr) 



Efe>0 c k 



(21) 



We have evaluated this sum numerically and plot- 
ted it against r, for a fixed value of £ in Fig. (4a), 
inset which shows its non-monotonic nature. In the 
continuum limit, the argument of arcsine takes the 
form [2 cos(2<) - 2vr 2 £ + 27r£Si(27r£) - 2twi7V(£, t)] 
where -^"(£, t) is defined as 



: exp 



(-JL) 



1 



which shows explicitly that GcD{r,t 1 L) is a function of 
^ and r only. 

To measure Gsp{r,t, L) for particles on an EW sur- 
face we performed Monte Carlo simulations as be- 
fore. After equilibrating the system, we measure 
-£ EiLi (J i(0)o'i +r (t) for about L z /10 time steps, then af- 
ter a gap of a few hundred time steps, we take another 
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FIG. 3: Aging auto-correlation for CD and SP models with 
(a) EW and (b) KPZ interfaces. The cusp exponent (3' was 
determined for ti ^> ti, after subtraction from m 2 . The CD 
model data in (a) has been multiplied by 1.5 to distinguish it 
from the SP model data points. The inset shows the power 
law behavior in the regime t\ < i2. We used L — 2048 in 
(a) and L = 8192 in (b). The Inset shows the data with 
h = 500,2000,8000 in both (a) and (b). For extraction of /3' , 
we used ti = 2000, 8000, 32000. 
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FIG. 4: 27ie time dependence of G(r, t, L) is shown for parti- 
cles on an (a)EW and (b) KPZ surface for £ = 0.016. The 
values of L are 256, 512, 1024 for (a) and 512, 1024, 2048 for 
(b). The scaled auto- correlation is also shown, for compar- 
ison. The insets show the same quantity calculated for the 
corresponding CD model with = 0.125 for both cases. 



set of data. We finally average over 10 5 such histories. 
The results are shown in Fig. (4a) where we have also in- 
cluded the scaling function h(r) to compare the long time 
behavior. The corresponding results for KPZ surface are 
shown in Fig. (4b). 

VII. LARGEST CLUSTER IN STEADY STATE 

One of the key characteristics of FDPO is the presence 
of strong fluctuations, even in the thermodynamic limit. 
In the steady state, large clusters are present in the sys- 
tem and the cluster size distribution follows a power law. 
As a result of fluctuations, these clusters undergo large 
changes in their lengths, associated with the fact that 
the macroscopic state of the system keeps changing. For 
a system of size L, the typical lifetime of a macrostate 
scales as L z . The question arises: if the lifetime of a state 
is so much smaller than exponential, in what sense can 
we call such a state a 'phase'? We have addressed this 



question in the following way. Let l m ax(t) be the length 
of the largest cluster present in the system at time t. In 
a disordered state, this length scales as logi. But start- 
ing from a random initial configuration, as the system 
approaches steady state, l ma x(t), although a fluctuating 
quantity, shows an increasing trend. Finally, in steady 
state, l m ax{t) is still fluctuating, thereby changing the 
macroscopic state of the system. But l ma x{t) continues 
to remain substantially above its disordered state value 
logL. In other words, the system manages to retain its 
ordered character despite steady state fluctuations. The 
system continues to move from one macroscopic state to 
other over a time-scale of L z . But each of these states 
are ordered in the sense that they all correspond to large 
values of l max (t). 

We have studied the distribution of l max in steady state 
as well as in disordered state. Our studies show that 
as system size increases, the overlap between these two 
distributions falls off. This means that as L grows, it is 
increasingly unlikely for the steady state l max to come 
down as low as its disordered state value. The time-scale 
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for such a transition in fact grows exponentially with L. 

The disordered state distribution was obtained by aver- 
aging over 10 8 data points. The mean of this distribution 
scales as log-L as mentioned earlier. 

After the system has reached steady state, we mea- 
sure the largest cluster present in that configuration. We 
let the configuration evolve in time and after waiting for 
few hundred time steps, we again measure l m ax{t)- We 
obtain the distribution P(l ma x, L) after normalizing over 
10 6 such data points. As shown in the following figures, 
P(l max , L) for different L values undergo a scaling col- 
lapse as l m ax is rescaled by the mean of the distribution 
(I). We have found that (I) ~ L^, where the exponent (f> 
depends on the dynamical rules. For particles on an EW 
surface <f> ~ 0.86, whereas for KPZ advection, (j> ~ 0.60 
while for KPZ anti-advection cj> ~ 0.91. We show the 
data for KPZ advection in fig.JSJ). 



0.06 




40 80 



■max 

FIG. 5: The distribution of the length of the largest cluster 
P(lmax, L) for particles advected by KPZ surface is shown for 
L = 256, 512, 1024, with the scaling collapse in the inset. The 
curves to the left show the same distribution in disordered 
state, after rescaling the y-axis by 0.2 



VIII. DISCUSSION 

In this paper, we have studied the dynamics of inter- 
acting passive scalars driven by a fluctuating Edwards- 
Wilkinson or Kardar-Parisi-Zhang surface (or equiva- 
lently, a Burgers fluid), by characterizing the scaling 
properties of correlation functions in steady state and 
those of aging correlations during the approach to steady 
state. It is instructive to compare our results with ear- 



lier work on the dynamics of passive scalars in different 
contexts. 

Mitra and Pandit 9] studied the dynamical properties 
of a system of non-interacting passive particles, advected 
by an incompressible fluid, whose velocity field is drawn 
from the Kraichnan ensemble, and therefore has power 
law correlations in space, but is delta-correlated in time. 
By contrast, we have studied passive particles with hard- 
core interactions, advected (in the Burgers case) by a 
compressible flow which has power law correlations in 
time. The compressible nature of the flow results in par- 
ticles being driven together in our case. The resulting 
state is described by a space-time correlation function 
G(r, t, L) which is a function of the scaling combinations 
r/L and t/L z , as in any phase-ordering system. However, 
there is no non-trivial r — t scaling in the limit L — > oo. 
On the contrary, in [9|, G(r,t,L) was found to show a 
scaling between r and t, with t ~ r z for fixed L. This 
difference of behavior is presumably a reflection of the 
strong differences between passive scalars with clustering 
or phase-ordering tendencies, and those which spread out 
in space. 

Even when the driving fluid is compressible, the de- 
gree of clustering of the passive particles depends on the 
nature of interactions between them. In the presence 
of hard-core interactions, the system reaches a phase- 
ordered state, albeit one with strong fluctuations. As a 
consequence, in the limit of small scaling argument, the 
spatial and temporal correlation functions show a cuspy 
approach to a finite intercept. However, in the absence 
of any interaction, the passive particles go into a more 
strongly clustered state, where the correlation functions 
show a power law divergence at the origin 5]. 

Finally, the study of the largest cluster allows us to 
arrive at a simple picture of a fluctuation dominated 
phase-ordered state. Despite the presence of strong fluc- 
tuations, the system never loses its ordered character. 
Fluctuations carry the system from one ordered config- 
uration to another macroscopically distinct one, over a 
time-scale ~ L z . However, the probability for the sys- 
tem to leave this attractor of ordered states vanishes ex- 
ponentially with the system size. 



IX. ACKNOWLEDGMENTS 

We acknowledge useful discussions with S.N. Majum- 
dar and S.Ramaswamy. SC would like to thank TIFR 
Endowment Fund for partial financial support. 



[1] B.I.Shraiman and E.D.Siggia Nature 405, 369 (2000). 
[2] E.Balkovsky, G.Falkovich and A.Fouxon Phys. Rev. Lett. 

86, 2790 (2000). 
[3] K.Gawedzki and M.Vergassola Physica D 138, 63 (2000). 
[4] E.Medina, T.Hwa, M.Kardar, Y.-C. Zhang Phys. Rev. A 



39, 3053 (1989). 
[5] A.Nagar, M.Barma and S.N.Majumdar Phys. Rev. Lett 

94, 240601 (2005). 
[6] B.Drossel and M.Kardar Phys. Rev. B 66, 195414 (2002) 
[7] D.Das, M.Barma and S.N.Majumdar Phys. Rev. E 64, 



9 



046126 (2001). 

[8] G.Manoj and M.Barm a J. Stat. Phys. 11 0, 1305 (2003). 
[9] D.Mitra and R.Pandit \nlin.CD/04120l3\ 
[10] S.F.Edwards, D.R.Wilkinson Proc. R. Soc. London A 
381, 17 (1982). 

[11] A.-L.Barabasi and H.E.Stanley, Fractal Concepts in 
Surface Growth (Cambridge University Press, New 
York,1995). 



[12] A.J.Bray Adv. Phys. 43, 357 (1994). 

[13] G.Porod, in Small Angle X-ray Scattering, edited by 

O.Glatter and L.Kratky (Academic Press, New York, 

1983). 

[14] P.Meakin, Fractals, Scaling and Growth Far from Equilib- 
rium (Cambridge University Press, Cambridge, 1998) 



